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Abstract. There is a fairly tight correlation between the pulsation periods and effective 
temperatures of ZZ Ceti stars (cooler stars have longer periods). This seems to fit the the- 
oretical picture, where driving occurs in the partial ionization zone, which lies deeper and 
deeper within the star as it cools. It is reasonable to assume that the pulsation periods should 
be related to the thermal timescale in the region where driving occurs. As that region sinks 
further down below the surface, that thermal timescale increases. Assuming this connection, 
the pulsation periods could provide an additional way to determine effective temperatures, 
independent of spectroscopy. We explore this idea and find that in practice, things are not 
so simple. 



1. Introduction 

ZZ Ceti stars, also called DAVs, are the 
coolest known pulsating wiiite dwarfs (see fig- 
ure 1). They show only hydrogen in their 
spectra and are observed to pulsate in a nar- 
row instability strip around 12000K. The num- 
ber of known ZZ Ceti stars recently dou- 
bled (Mukadam 2004), and more have been 
found since (Mullally et al. 2005, Kepler et 
al. 2005, Castanheira et al. 2006, Castanheira 
2006). With such a significant number of stars 
known, it was promising to study their ensem- 
ble characteristics in more depth (Mukadam 
2004, Mukadam et al. 2005). 

Even before the new ZZ Ceti's were dis- 
covered, a clear relation between their mean 
pulsating periods and spectroscopically deter- 
mined effective temperatures was already ap- 
parent (Clemens 1993, Kanaan 1996). Cooler 
ZZ Ceti's show longer periods (see figure 3). 
Well aware of the dificulties associated with 



the spectroscopy of such stars, Mukadam et al. 
(2005) suggested that we could use this rela- 
tion to obtain an independent measurement of 
temperatures for those stars. 

We compared relevant thermal timescales 
of models with different effective temperatures. 
We looked in particular at driving through the 
K-y mechanism (Cox & Giuli 1968) and con- 
vective driving (Brickill 1991, Goldreich and 
Wu 1999). We expected that if we picked the 
right convective efficiency, the mean periods 
would agree with the thermal timescales (to 
within a multiplicative constant). We would 
then use the effective temperatures of the mod- 
els as a temperature scale. 

Besides helping the ongoing effort to bet- 
ter define the ZZ Ceti instability strip, we then 
planned to apply the same method to the less 
numerous DBV stars (Helium atmospheres). 
Spectroscopic determinations of effective tem- 
peratures for those stars are difficult to do and 
uncertainties are large. We need more accurate 
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temperature determinations for DBV stars be- 
cause they can be used to learn about plasmon 
neutrinos, and the results are most sensitive to 
temperatures (Winget et al. 2004, Kim et al. 
2005, Kim et al. 2006) 
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Fig. 1. H-R diagram showing the different classes 
of pulsating stars. In this article, we are interested in 
the DAV's. 



2. Pulsation driving and timescales 

Pulsations in ZZ Ceti stars are driven in the 
partial hydrogen ionization zone. We looked 
at two processes, each of which are associated 
with a thermal timescale. The first process is 
the K-y mechanism, and the relevant timescale 
in this case is the thermal timescale at the base 
of the partial ionization zone {Rb). Like Winget 
et al. (1982), we assume that only modes that 
have periods of the same order as the thermal 
timescale are driven. This thermal timescale is 
the time it takes to transport energy from Ri, to 
the surface of the stai\ It is equal to the heat 



content of the layers above Ri, divided by the 
total luminosity: 

^MrCvT 

Tth (1) 

The second process is convective driving. 
Convective driving was proposed to attempt 
to resolve a discrepancy between the peri- 
ods of the overstable modes and the above 
timescale. Noting that the convective turnover 
time was much shorter than the periods of the 
modes, Brickill (1991) assumed that convec- 
tion responded instantanuously to the pulsa- 
tions. Goldreich and Wu (1999) later revisited 
his theory of convective driving and arrived 
at a new thermal timescale, which is the time 
during which the convection zone can bottle 
up the flow that enters from below. This new 
timescale Tc diff'ers from the traditional thermal 
timescale by a multiplicative constant roughly 
equal to 4. 

So we have two timescales for which we 
know the dependence on the local tempera- 
ture. To compare with the observed periods, we 
need them in terms of effective temperature. 
Numerical calculations (Montgomery 2005) 
show that Tc, and therefore also t,/, are propor- 
tional to the effective temperature to a power 
roughly equal to -90 (for DAVs). This re- 
sult is insensitive to the convective efficiency. 
Therefore in a log P - log T^b plane, we would 
expect the data to lie on a line with a slope 
of -90, no matter which of the two theories is 
the correct theory for driving in ZZ Ceti stars, 
and independently of the right convective effi- 
ciency. 

3. Results 

Figure 2 shows the data on a log P - log T^b 
plane along with curves that show rdTf^ff) for 
different convective efficiencies. We see that 
indeed, the thermal timescale curves are par- 
allel to each other, and we also see that they 
have a different slope than the data. 

Figure 3 shows a subset of the data with er- 
ror bars. We have retained stars in a narrower 
\ogg range (8.0 to 8.2) in a rudimentary at- 
tempt to remove the logg dependence of the 
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Fig. 2. The ZZ Ceti stars discovered before 
Mukadam et al. 2004 (filled circles), and the ones 
found by Mukadam et al. from SDSS data (open cir- 
cles), along with for different convective efficien- 
cies. In the legend, we designate Bohm-Vitense's 
treatment of convection (Bohm-Vitense 1958) by 
MLl, and the more efficient Bohm & Cassinelli con- 
vection (Bohm & Cassinelli 1971) by ML2. a is 
the ratio of the mixing length to the pressure scale 
height. There is a systematic offset of ~ 300K be- 
tween the two samples, probably caused by system- 
atic differences between Bergeron's models atmo- 
spheres and Koester's, for which we artificially cor- 
rected by substracting 294 K from the "old" sample. 
The choice is somewhat arbitrary since we are only 
concerned with relative temperatures at this point. 
The three open circles at high effective temperature 
and high < P > that seem to deviate from the cor- 
relation also happen to have high gravities. This is 
consistent with theory, which predicts that r^. is sig- 
nificantly higher for high gravities (about an order 
of magnitude for Alog(g) = 0.6) 



timescales. We fit one line to the cool DAVs, 
one to the hot DAVs, and a third to the full ef- 
fective temperature range. The boundary is at 
11320K. The slopes of the three lines are, re- 
spectively, -31, -61, and -28. Uncertainties are 
large and those results are not entirely incom- 
patible with the expected slope of -90. The dif- 
ference in slope between the red-edge pulsators 
and the blue-edge pulsators is not statistically 
significant (again because of the uncertainties) 
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Fig. 3. A subset of the ZZ Ceti stars (those with 
logg between 8.0 and 8.2) and different fits. The 
black line is a fit to the whole subset (slope = -28), 
the thin grey (orange) line is a fit to the cool end stars 
(slope = -31), and the thick grey (light blue) line is 
a fit to the hot end stars (slope = -61). The boundary 
between the cool (cDAVs) and hot stars (hDAVs) is 
11320 K. 

but theoretically the blue-edge pulsators don't 
show significant convection zone. 

We do note, however, that departures from 
the models become important for cooler stars. 
It looks like a saturation mechanism is at work, 
which prevents modes much higher than 1000s 
from being driven, though that is not the only 
explanation. 

4. Some ideas 

In our analysis, we have made a few assump- 
tions. First we assumed that the convective ef- 
ficiency remains the same as the stars cools. 
Maybe convection becomes less effective with 
decreasing effective temperature, which would 
explain the observed shorter timescales. 

While this may seem to go against com- 
mon wisdom, there are a few physical pro- 
cesses which could work against convection 
for cooler stars. Magnetic fields may become 
more important as the star cools and the depth 
of the convection zone increases. One could 
easily imagine the magnetic field lines con- 
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straining the motion of the plasma and inhibit- 
ing the convection, leading to a lower convec- 
tive efficiency. 

Most red edge ZZ Ceti stars show large am- 
plitude pulsations. The linear approximations 
break down, and one could easily imagine that 
such high amplitude pulsations could have an 
effect on the convective efficiency. At the red 
edge, the pulsation amplitudes can be as high 
as 20% at the surface (Kleinman 1995). At 
those amplitudes, the energy carried by the pul- 
sations is non-negligible, which would lead to 
a less efficient convection. This of course fails 
to explain why the red edge ZZ Ceti's that are 
not high amplitude pulsators deviate from the 
theoretical thermal timescales in the same way 
the others do. 

Maybe it is our assumption that the 
weighted mean period (i.e. the period(s) of 
the dominant mode(s)) is ~ Tc that is flawed. 
Again, this may not work so well where 
non-linearities are important. At lower effec- 
tive temperature, the convection zone grows 
deeper, and when deep enough, maybe its mo- 
tion affects the timescales. 

The most obvious explanation, and maybe 
the best one we have is that modes with peri- 
ods higher than about 1200 seconds are sup- 
pressed through a different mechanism. Lack 
of surface reflection comes to mind (Hansen et 
al. 1985). Modes with long wavelengths (long 
periods) are not reflected very efficiently in the 
outer layers of the stars where the density falls 
off exponentially, and lose their energy. This 
places an upper limit on the periods of the 
driven modes. Hansen et al. (1985) did a sim- 
ple non-adiabatic calculation and found that 
modes with periods greater than 3000s could 
not sustain themselves in ZZ Cetis, although it 
is possible that The actual limit may be more 
Hke 1200s. 

5. Conclusion 

We set out to develop an independent method 
to determine temperatures for ZZ Ceti stars, 
with the intention to apply the same method 
later to DBV stars. We based our analysis on 
the assumtpion that the period(s) of the dom- 
inant mode(s) should be at least proportional 



to timescales corresponding to several theories 
of driving in ZZ Ceti stars. In actually doing 
the calculations, we found that this was not 
the case. While the theory predicts that in the 
log < P > - logTeff plane, the data should lie 

on a line with slope 90. We find the actual 

slope to be less steep. 

We give a few ideas to explain the discrep- 
ancy, but this remains somewhat of a mystery, 
which would be very interesting to elucidate in 
order to not only obtain better effective temper- 
atures, but also to improve our understanding 
of driving in ZZ Ceti stars. 
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